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1 Introduction 


Fermi ball(F-ball), which is a kind of nontopological solitons, is hrst proposed as a can¬ 
didate for cold dark matter(CDM) |]I|. Similar object is considered to have possibility to 
explain the baryon nnmber asymmetry of the present nniverse |^. An F-ball consists of 
a closed domain wall and fermions which localize on the wall. Such an object is thought 
to be produced after a phase transition, in which two almost degenerate vacua exist. At 
this phase transition, two regions, which correspond to two vacua, and domain walls are 
produced. If these two vacua are completely degenerate, the domain walls will dominate 
the energy density of the universe soon |Q]. However, if the energy density of one vacuum 
is slightly smaller than the other one, the true vacuum pushes the domain walls and re¬ 
duces the false vacuum region. If some fermions are captured on the domain walls, the 
left region may be stabilized due to the Fermi pressure of the captured fermions. 

We hrst assume that the region has complete spherical form with radius, R. In this 
case, the total energy of this system consists of three parts, the surface, the volume, and 
the Fermi energy within the thin wall approximation: 

= E, + E, + Ef = ( 1 ) 

O Oit 

Here, S, e, and N are the surface tension, the energy density difference between the two 
vacua, and the number of fermions on the wall, respectively. Since the surface energy and 
the volume energy are the increasing functions of the radius and the Fermi energy is the 
decreasing one, this object is stabilized at a certain radius. We deal the volume energy 
as a perturbation term hereafter, since it is much smaller than the other energies in most 
cases. Minimizing this total energy with respect to R, we hnd out that it is proportional 
to the number of the fermions on the wall with the following critical radius, Rc- 

Etotai = K = (127rS)3). (2) 

d K 

This object is called as Fermi ball. 

If this large F-ball were stableQ, it were interesting as a candidate for CDM. It, however, 
is unstable against deformation from spherical shape and fragments into small pieces |T[]. 
This can be understood as follows. Since the energy of an F-ball is proportional to the 
number of the fermions on the wall as above, it has same energy even if it is divided into 
some smaller F-balls. However, there exists the small volume energy, which we neglected. 
Taking into account this contribution, we can easily hnd out that the fragmented state has 
smaller energy. An F-ball therefore continues to be small until the thickness of the wall 
becomes comparable to its radius and the thin wall approximation breaks down. Such 
an F-ball is too small and there must be large number of them in order that they have 

* We call an F-ball large if its radius is much larger than the thickness of the wall. 
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sizable contribution to CDM. An F-ball was, therefore, thought to consist of a domain 
wall and a fermion which interact with ordinary matter only very weakly [Q. 

An electrically charged F-ball, which is stable against the fragmentation owing to the 
long range electric force, is then proposed [Q]. This F-ball can be large and heavy enough 
to have sizable contribution to CDM even if they are very dilute. 

We, however, have three questions about the above treatments of an F-ball: 

• Is an electrically neutral F-ball unstable even if the curvature effect is taking into 
account? 


• Is an electrically charged F-ball really stable? 

• Is an electrically charged F-ball stable even in the hot early universe where the 
F-ball is consider to be produced? 


First, when we consider an electrically neutral F-ball, we used the thin wall approx¬ 
imation and neglected the contribution from the curvature of the wall. Though this 
contribution is really much smaller than the total energy of an F-ball, it can be very 
important. Neglecting the curvature effect, we concluded that the energy of an F-ball is 
same as that of its fragmented state in the absence of the volume energy. However, if 
we take into account the contribution from the curvature, two states will have different 
energies. This curvature contribution, therefore, would determine which state is stable. 

Second, The total energy of an electrically charged F-ball approximately consists of 
the surface energy and the coulomb energy for large N because the coulomb energy is 
proportional to N'^: 


Etotal — T A'c 




aN'^ 


(3) 


Minimizing the total energy with respect to the radius, we obtain the total energy as 
follows: 


^total 


1 

3 




(4) 


4 

Since the total energy is proportional to Ns, a fragmented state has apparently smaller 
energy than an F-ball. We therefore can not conclude that the electric force stabilize the 
F-ball from the fragmentation only by comparing the energies of the two states. We then 
should investigate that is there any energy barrier to deform an F-ball to a fragmented 
state. If there is enough barrier due to the electric long range force, an F-ball can be 
metastable. 

Third, in thermal bath, the electric long range force becomes short range one with 
typical Debye screening length. The energy barrier, which arises from the electric long 
range force, would become week or disappear due to this Debye screening. 

To answer these three questions is the main aim of the present paper. Details of the 
following analysis are shown in Ref.[Q. 
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2 Stability of a neutral F-ball 


We first consider the curvature effect to the stability of an F-ball. Though we only consider 
a simple model described by the following Lagrangian density, most results will be able 
to be applied to more complicated models [^: 

c = liMf + Mira, - G4 >»f - U{4>) , (5) 

where 0 and 'ipp are a scalar and a fermion held, respectively, and G is a Yukawa coupling 
constant. Here, t/(0) is an approximate double-well potential i, 

UW = ^ - «")" + U.W . ( 6 ) 

taking the second term, which breaks Z 2 symmetry under 0 — 0 , much smaller than the 

hrst one. This model has a kink solution^], which interpolates the two vacua, 4> = v, —v: 

(l){z) = V tanh ^ = 2/(A/Ar;)). (7) 

On 

We assume the width of the wall, 6n is much smaller than the F-ball size. We call the 
plane where 0 vanishes ’’surface” of the domain wall in the present paper. This domain 
wall has the surface tension, S: 


S = 


271 o 

■ 


( 8 ) 


Substituting Eq.(|^ into the Euler-Lagrange equation of 'ipp, we hud that 0^? has a zero¬ 
mode solution, 

0^(0 = 0ir(2;o)exp|-G^ dY 0(Y)| . (9) 

Here, 0i?(zo) is a spinor eigenstate of *73 = —1. A fermion can localize on the domain 
wall owing to this zero mode. This model is therefore a suitable model to investigate the 
properties of an F-ball, since it has the domain wall solution and the fermion, which can 
localize on the wall. 

We evaluate the next-to-leading order approximation of the thin wall expansion in 
this model. We assume that the scalar held expectation value does not change along the 
surface and write the surface energy as follows^: 

E00] = yd^x |^(V0(x))2 + ^ (0(x)2 

^Note that the qualitative discussions in the following do not depend on the explicit form of U{(/)). 

^ Though we neglect the effect of Ue to show the following explicit solution, the presence of the solution 
is not affected by Ue- 

^ We only consider cases, i?i, i ?2 > 0, since it is sufficient to analyze the stability of an F-ball. 
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where 0 stands for the derivative of 0 with respect to w. The local coordinate, w is 
perpendicnlar to the snrface and becomes zero on the snrface. The cnrvatnre radinses, 
Ri and R 2 are principle cnrvatnre radinses on each points and the local coordinates, 61 
and 62 are the corresponding angnlar coordinates (see Fig.^. 



Fignre 1: Local coordinates on the wall. 
Taking variation of with respect to 0, we obtain the equation, 

0 + ( p ^ + T.\ ) 0 = ^0(0^ - 

\Ri + w R 2 + WJ 2 
with the boundary condition. 


( 12 ) 


±v {w —>• ±CX)) 

0 (w = 0) . 


(13) 


We now estimate Eg within the thin-wall approximation, ignoring the curvature effect.]^ 


The equation (|T^) becomes in the leading order as. 


^00 (00 




(14) 


with the boundary condition. 


±v {w —>• ±cxd) 

0 (w = 0) . 

^This of course leads to the same result as Macpherson and Campbell derived. Our derivation is how¬ 
ever meaningful as a step toward the estimation of the energy up to the next-to-leading order contribution 
in the thin-wall expansion. 
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The solution to Eq.(0) is the kink, 


... , tc 

(Po[w) = utanh — . 

On 

(16) 

Adding the Fermi energy, we obtain the total energy of the neutral F-ball within the 
thin wall approximation, E^: 

E% = J:S + ^JdS , 

(17) 

where S is the surface tension in Fq.(§). We should minimize 
the total fermion number. 

this total energy, keeping 

N = J dS npi^p) , 

(18) 

constant. Using the Lagrange’s multiplier method, we hnd that npi^p) is constant: 

ufM = • 

(19) 

From this, we hnally get the total energy: 


Sv’S ■ 

(20) 

Minimizing with respect to S, we obtain. 


E% = (127rS)^/3jY ^ 

(21) 

with the critical area of the surface. 



(22) 


Since the total energy is proportional to the number of the fermion, we can not know the 
stability of the F-ball in this order approximation. 

We then calculate the next-to-leading order approximation of the thin wall expansion. 
We here expand 0 and E]^ with respect to Sjsi/Ri and Sn/R 2 , 

0 = 00 + 01 H- 

E]\f = E^ E^ E"^ -!-■■■ . (23) 

Since E]^ vanishes owing to the variational principal, we estimate the energy up to Ej^. 
It is enough to expand 0 up to 0i in this case. From Fq.(|T^), 0i satishes, 

01 - ^ (300 - 01 = “ 00 ’ 
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with the boundary condition, 


01 


0 (ic —>• ±cx)) 
0 (m; = 0) . 


The solution to Eq.(^) is, 




1 

2\/X V-Ri R 2 ) 


Mw/Sn) 


1 


with. 


fi{x) = sech^x — cosh^ x + -sinh^xtanh^x 

o 


sech^x 


|sinh4x + 8sinh2x + 12x| . 


Substituting 0o and 0i into Eat, we obtain. 


El = Cl, dS 




with. 


^2 _ 

Wat — 


1 

7X. 


Ch — X— j dy, ( 


4V^ 


w‘ 4 >o~ - '^ 0 / 1 ) *' 

' . ; . 0.28!) 
dTO Wi ^-^ < 0 . 


(25) 


(26) 


(27) 


(28) 


(29) 


We consider the stability against the fragmentation from the next-to-leading order con¬ 
tribution, E^. The hrst term in Eg . (12^) is proportional to the integration of Gaussian 
curvature on the closed surface, which is known to be dvr (Gauss-Bonnet’s Theorem). We 
therefore hnd out that the hrst term does not depend on the shape of the F-ball and 
depend only on the number of the F-ball. Since is negative, a fragmented state has 
smaller energy. On the other hand, the second term is zero only for sphere and negative 
for any other shapes. A spherical F-ball then deforms and fragments to small pieces. We 
conclude that a neutral F-ball is unstable even in the absence of the volume energy due 
to the curvature effect. 


3 Stability of an electrically charged F-ball 

We consider the stability of an electrically charged F-ball in thermal bath in the present 
section. We here consider the free energy of the fermion gas on the wall and the electron 
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gas, which surrounds the F-ball|: 

Fc\v,] = /d’x |-i(vv;(x))V^,(v;(x)) 

- nF(xF)U(x)d(|x - xf|)| , (30) 

where jFg is the free energy density of the electron gas. When the temperature, T is much 
higher than the electron mass, Te can be written as. 


-^e(K) = 


2T 


(2vr)' 


r^ 2 y 2 

* e 

6 


d’p I f + log (l + 


T 


P-Ve 

e ^ 


^ e 

127r2 


+ log(l + e 

+ (Ve independent terms). 


(31) 


The free energy can be obtained by extremizing Eg . (1301) with respect to 14 (x). Since 
14 (x) is expected to change rapidly near the wall, we also use the thin wall expansion 
here. Assuming that 14 depends only on tc, we can express Fq as. 


Ff| 14| = /d5d»(l + |^)(l + f){-^K 


V 2 

e ^ 


^ - nFV,S(w) 


(32) 


Taking variation of Eq.(|3^ with respect to w, we obtain the equation. 




+ 


e"^ \Ri + w R 2 + w 




TWe ^ W 


Svr^ 


+ np5{w) 


with the boundary condition. 


(33) 


V,^0 


[w 


±CX)) 


(34) 


This corresponds to the well-known Thomas-Fermi equation. 

We hrst calculate Fq within the thin-wall approximation, ignoring the curvature effect. 
From Eq.(^), we obtain the equation. 


- 1/0 

=2 e 


T2I4O 

3 


+ 


{Ky 

37r2 


npdiw) , 


(35) 


with the boundary condition. 




0 


[w 


±CX)) . 


(36) 


II Though we assume that the fermion on the wall has an electric charge, we can easily extent our 
analysis to other cases. 
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The solution of Eq.( 


with 


IS 


V}>(w) = 


-V2ttT 


sinh {Q + ct 
Vs 


At — 


eT ’ 


and with, 


ct = cosh 


.1 1 + vT+^ 


a = 


a 


VSerip 

V27iT^ 


(37) 


(38) 


(39) 


Substituting Eq.(p7D into Fc, we obtain, 




2V27ie 

\ sVs Vs ■ 


(40) 


Adding the surface energy, we obtain the free energy of the F-fall within this order ap¬ 
proximation. 


= S5 + 


iy/n ^ 2A/27re\ 

3 v^ / VM ' 


+ 


(41) 


Since the coulomb energy only change the coefficient of the second term, we can not know 
the stability of the F-ball within this order approximation, too. We therefore need to 
calculate the next-to-leading order approximation. 

We expand W and Fq with respect to Sc/Ri and Sc/R 2 0, 


n = n“ + K' + --' 

Fc = i^ + F^ + Fa + '" ■ 

Since Fq vanishes, we estimate F^, expanding W up to VV From Eq. 


(42) 


we obtain 


1 

—W - — + 


(vSf 


TT^ 


=-( — + — ) FO 

e _9 I D D < e ’ 

rxi 1X2 


with the boundary condition. 


—> 0 (w —>• ±cxd) 


(43) 


( 44 ) 


*Here, 5c is a typical screening length near the surface of the F-ball, 5c — XtCt- 














The solution to Eq.(^) is, 


Here, / 2 (|w|/Ar) is defined as, 

fn \/\ \ w y/Qir cosh {\w\/Xt + ct) r, ,,, , . f ^ 

/2(f|/At) - --p-r-—T7]—rTT—^-r I/sIfI/At + ct) - /sIct)} 

\w\ e smh {\w\/Xt + ct) 

where fsix) is the following function: 

1 2 1 

hix) = -sinhxcosha: H—tanhx-sinh^x — x . 

jsv y 3 3 3 

Substituting and into Fc, we obtain, 

1 


F^ = Ch I dS 


RiR 2^^^ J \Rl R2 


with. 


ch = 


(^2 _ _ 

“ 4e2 


r J dw - ^0/2} ~ - 27 r "/3 ^1/6^ < g 

r . / 9 \ 

J dw V°f2 ~ (^) ^ 


Adding the surface energy, we finally obtain. 


F^ = Ei + F^ 


= dS 


—— + C^ dS — - — 

R1R2 J \Ri R2 


( 45 ) 


(46) 


(47) 


(48) 


(49) 


(50) 


where C* (i = 1, 2) is the sum of and C^. Since the first term is always negative and 
depends only on the number of the F-ball, a fragmented state has smaller energy than an 
F-ball state. However, the second term is positive in range > 0, where is. 


A^e 4 . 


(51) 


In this range, an F-ball is stable against the deformation from spherical shape. It, however, 
will fragment due to the tunneling effect or the thermal fluctuation, since a fragmented 
state has smaller energy than an F-ballQ. We therefore conclude that an electrically 
charged F-ball is metastable in this range. 

t^The condition that an F-ball survives until now will be shown in Ref. [Q. 
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4 Summary and Discussion 


We considered the stability of an F-ball. We found out that an electrically neutral F-ball 
is unstable even if we take into account the curvature effect. We notice that this is true 
even in the absence of the volume energy. We then found out that an electrically charged 
F-ball is metastable even in thermal bath for the parameter range, Eq.([^). 

We lastly discuss the allowed region of parameters, e, N, and n. For the volume energy, 
which arises from the energy density difference between the two almost degenerate vacua, 
e, we have two constraints from cosmology and from the stability of the electrically charged 
F-ball. From cosmology, the energy density difference should be large enough to avoid 
the black hole dominated universe 0: 


This condition is rewritten as. 





(52) 


(53) 


On the other hand, the energy barrier between an F-ball state and a fragmented state 
should not vanish due to the volume energy for the metastability of a charged F-ball. In 
order to satisfy this condition, we have. 




iTlR^ 


This condition is rewritten as. 


3a 2 >( 


(54) 


(55) 


We next consider a constraint from experiments. Assuming that F-balls have sizable 
contribution to CDM, we obtained a constraint for the energy of an F-ball in Ref. ^ as, 

E>10‘^^[GeV] (fi;>10[GeR]). (56) 


This condition is rewritten as. 


N> 


1025 - 

K[GeV]_ 


(57) 


We hnally consider how large the symmetry breaking scale, k is? In order that the 
conditions, Eq.(|53D and Eq. (|55[) , are compatible, we have. 


1 

■3a’ 


Iat- 


l> 



(58) 
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This condition is rewritten as, 



(59) 


In order that the conditions, Eg .(^71) and Eq.(^), are compatible, we have, 

[GeE]. 


(60) 


The symmetry breaking scale seems not to be so far from the electroweak scale. These 
constraints would help us to make a realistic model of an F-ball. 
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